The complex sine-Gordon model on a half line 
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In this paper, we examine the complex sine-Gordon model in the presence of a 
boundary, and derive boundary conditions that preserve integrability. We present 
soliton and breather solutions, investigate the scattering of particles and solitons 
off the boundary and examine the existence of classical solutions corresponding to 
boundary bound states. 



1 Introduction 

Two dimensional integrable field theories have become an area of extensive study. Their 
rich underlying mathematical structure allows for their exact solution. This in turn pro- 
vides valuable information about the wide range of physical phenomena which integrable 
field theories can be used to model, and more generally about non-perturbative field the- 
ory. A number of models have been studied in the presence of a boundary. This has led 
to results for both the classical and quantum scattering of objects, like particles and soli- 
tons, off the boundary. In particular Toda models |lj, the sine-Gordon model 0, and the 
sinh-Gordon model amongst others have been studied on the half-line. 

In this paper we study the classical two dimensional complex sine-Gordon (CSG) model 
in the presence of a boundary. The model in the bulk is described by the following La- 
grangian 

. ldudu* + dudu* „ , 

Ccsg = 77 ; ; , 1 

2 1 — £uu* 
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where d = d/dz, d = d/dz, z = (t — x)/2, z = (t + x)/2 . The field u is complex, (3 is a real 
coupling constant and the parameter £ can be changed by rescaling u and will henceforth 
be set to one. The Lagrangian is invariant under global phase rotations of the field u and 
this leads to a conserved U(l) charge. 

As such, the complex sine-Gordon model comprises an integrable generalisation of the 
sine-Gordon theory, with an internal U(l) degree of freedom. It was first introduced 
independently by Lund and Regge as a model of relativistic vortices in a superfluid || ^| 
and by Pohlmeyer in a dimensional reduction of a 0(4) non-linear cx-model [Q. It belongs 
to a class of generalisations of the sine-Gordon theory, which in the literature is referred 
to as homogeneous sine-Gordon theories. The latter appear as G/U(l) gauged Wess- 
Zumino-Witten models perturbed by a potential. These generalisations describe integrable 
perturbations of c > 1 conformal field theories and have been studied in the bulk in j7|, Q . 
The CSG model appears as the simplest case where G = 577(2) and describes integrable 
perturbations of parafermions by the first thermal operator []T2l [L9|]. The quantum case 
was also studied by Dorey and Hollowood in fl8| , and Maillet and de Vega in [16 . 



The CSG model has found many applications in different fields of physics, from general 



relativity M, to the description of propagating optical pulses in a non-linear medium [10 



In the latter, the CSG theory was used as a generalisation of the pioneering work of McCall 
and Hahn [11], where the simple sine-Gordon (SG) was used as a field theory description. 
The CSG theory provides for a more realistic description of optical pulses incorporating 
effects like frequency detuning and modulation, while at the same time the theory can be 
extended beyond the description of two-level atom systems to multi-level atom systems. 

In the next section of this paper some aspects of the theory in the bulk will be reviewed. 
Classical solutions, including particles, solitons and breathers will be presented for the 
model in the bulk clarifying previous treatments. Unlike the sine-Gordon model where the 
soliton solutions are topological in nature, the CSG model solitons which carry a Noether 
U(l) charge, have no obvious topological charge associated with them. This follows from 
the trivial vacuum structure of the model for (3 > 0, although as we shall explain chargeless 
solitons do carry the topological charge of the sine-Gordon model in a subtle way. We also 
discuss the existence of breathers and find a class of charged breather solutions. 

In section 3, we consider the effect of introducing a boundary and suitable boundary 
conditions are found which preserve the integrability of the model. In order to do this we 
explicitly construct low-spin conserved charges by abelianising the Lax pair description 
of the model, and ensure that these are conserved on the half-line by imposing boundary 
conditions. We find a natural generalisation of the boundary conditions that Ghoshal and 
Zamolodchikov PUJ introduced for the sine-Gordon theory. 



The final section in the main body of the paper deals with the scattering of particles 
and solitons. We find the corresponding solutions and use them to calculate the classical 
time-delay. We also find solutions corresponding to boundary bound states. The paper 
concludes with some general remarks about the results presented, as well as a few open 
questions that are related with some interesting properties of the model. 
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2 Classical aspects of the CSG theory 



In this section we shall review some important features of the CSG theory and establish 
our notation. We present the equation of motion, the relation of the model with the sine- 
Gordon theory and review the mathematical background in which the CSG model can be 
regarded as a perturbed gauged Wess-Zumino-Witten model. Moreover vacuum, soliton 
and multi-soliton solutions are written down in a compact form. Finally, we clarify the 
existing confusion concerning the existence of breathers within this model and present 
explicit breather solutions. 



2.1 Definition of the model 

The Lagrangian of the model was presented in ([]]). We assume without loss of generality 
£ = 1 since the parameter can be dropped by rescaling the field variables. The equation 
of motion that follows is 

dBu + ^duBu + x _ =Q 
1 — uu* 

The relation between this model and sine-Gordon becomes obvious if we substitute 

u = smct)e 2ir) , (3) 

with (p and ry real fields. We use this notation here which gives a more standard connection 
with the sine-Gordon theory but differs from previous treatments by <fi — > | — 0. This 
takes ([!]) to 

£ = d<pd(p + 4 tan 2 (pdrjdr] - 4/5 sin 2 . (4) 

This is essentially the form of the Lagrangian derived by Lund and Regge, and Pohlmeyer. 
By taking the field rj to be constant the sine-Gordon Lagrangian emerges. As was demon- 



strated by Bakas [12] the theory can be reformulated in terms of a gauged WZW action. 



The corresponding action principle can be written as 

S = Sgwzw + Spot ■ (5) 
The action term S g wzw is the well known gauged WZW action 

Swzw = -— [ dzdzTr(g~ 1 dgg'' 1 dg) - — [ Ti^g^dg A g~ l dg A g~ x dg) 
47T Js LZir Jb 

+ — / Tri-Wdgg- 1 + Wg~ x dg + WgWg' 1 - WW) . (6) 

This action is defined in a three-dimensional manifold B whose boundary is our compact- 
ified normal two-dimensional space E. The field g is an 577(2) group element and g is the 
extension of g to the three dimensional manifold. The last term introduces gauge fields W 
and W which act as Lagrange multipliers. The S pot term is 

Spot = 7T f ^(g^sg" 1 ^) ■ (7) 
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This term breaks conformal invariance and thus gives rise to massive states. Varying the 
action yields the CSG equations of motion which can be expressed in a zero curvature form 



[d + {g- x dg + g-'Wg + i(3Xa 3 ) , 3 + (W - jg'^g) } = . (8) 
From the variation of the gauge fields W and W, two constraint equations arise 

dgg- 1 - gWg- 1 + W = , 

(9) 

g- 1 dg + g- 1 Wg-W = 0, 

which are critical in order to make the identification with the CSG theory. The connexion 
between the SU(2) matrix g and the complex field u of @ is given by 

^(-*r:*)< (io » 

where v = — uu*e~ ld . The field variable 9 should not be considered as an independent 
field but rather as an auxiliary field that is properly defined up to a constant through the 
constraint equations @. In the gauge where W = W = , the constraint equations take 
the form 

,u*du — udu* ^ ,udu* — u*du 

Ou = — % — ; — , ad = — % — — , 11) 

2(1 -uu*) 2(1 -uu*) v ; 

whilst the equation of motion now becomes 

[d-A,d-A] = 0, (12) 

where 

A = -(g- 1 dg + ip\<T 3 ) , A= l -g- l a z g . (13) 

This compact zero-curvature form of the equations of motion demonstrates the inte- 
grability of the model and will prove useful when we come to consider Backlund transfor- 
mations and conserved quantities in later sections. 



2.2 Vacuum solutions in the bulk. 

From ([]]) it is easy to see that the energy of the CSG model in the bulk is 

r ( \d u\ 2 + \d lU \ 2 A 

T-Cbuik = dx\ VAfJuu . (14) 

J \ 1 — uu J 

The most suitable candidate for a vacuum, would be a constant value for the field u that 
would force the kinetic term involving derivatives to vanish and at the same time minimize 
the potential term. It is clear to see from (^) and (|]) that for (3 > the obvious choice is 
u = 0, while for (3 < the choice should be \u\ = 1 if we insist that \u\ < 1 . The sign 
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of the parameter ft divides the theory into two sectors. In the matrix potential formalism, 
both sectors are treated simultaneously as the diagonal and off-diagonal parts of the field 
variable g. In this context the fields u and v are both solutions to the CSG equation each 
derived for a specific choice of ft, and each corresponding to a different vacuum . This is 
because the two sectors are connected by a duality transform which interchanges the sign 
of the coupling constant ft [[14j] and simultaneously interchanges the role of u and v . That 



is to say the theory is invariant under the change 

the latter representing a transform akin to the Krammers-Wannier duality of the Z n 



parafermion theory ||15|| . Taking into account the invariance of the theory under this du- 
ality transform, we shall concentrate in this paper on the j3 > sector which corresponds 
to the diagonal part of the matrix formalism. A suitable vacuum solution would be 

9vac = ( ie m \ ) • (16) 

This selection is consistent with the choices appearing in the beginning of this section with 
the diagonal ft > sector, corresponding to the u = vacuum , while the off-diagonal 
ft < to \v\ = 1. It is noted that the apparent singular behaviour of the Lagrangian at 
\u\ = 1 does not appear as a problem embedded in the theory but is a direct consequence 
of the fact that the gauge fields W, W are ill defined at the specific point. 

2.3 Spectrum of the model 

The CSG model, like the sine-Gordon theory, possesses both particle and soliton solutions. 
When small perturbations around the vacuum are considered 

u = + e(x,t) , (17) 

the theory becomes linear when higher order terms in e are ignored 

(8% -df)e(x,t) +m 2 e(x,t) = 0, (18) 

where m 2 = Aft. The solution to the above equation is the familiar plane waves solution 

e(x, t) = e- iwt (Ae ikx + Be~ ikx ) , (19) 

where k and uj are related through 

uj 2 = k 2 + m 2 . (20) 

Different techniques have been used for the construction of soliton solutions like the 
inverse scattering method [IB| and the Hirota method [O. However both methods yield 



results that are both cumbersome and difficult to manipulate. The Backlund transforma- 
tion for the CSG model provides a more elegant way to obtain soliton solutions and can 
be written in terms of two matrix variables g and / |14 
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g-^dg-f-'df-^Llg-'asf,^} = 0, (21) 



dgg-^-a.dfr' + ^igf-^-a.gf- 1 ) = 0. (22) 

It is easy to show that both / and g satisfy the CSG equation as well as the constraint 
equation in the specific gauge choice. Taking / to be an already known solution, one can 
generate a new solution through the equations presented above. One-soliton solutions can 
be derived by applying the Backhand transformation on the vacuum solutions g vac of (|T6|). 
Each sector of the theory provides us with two sets of two first order differential equations 
that can be integrated, in order to provide the one-soliton solutions. The diagonal elements 
of g which correspond to the /3 > sector give 



uu* = 



(23) 



t-Vx 



The one-soliton solution that emerges is 

cos(a) exp (2i\f]3 sin 

cosh (2 a//? ^o V u, j y / 1 _ v z 

where V and a are real parameters associated with the velocity and charge of the soliton 
respectively. This solution was originally derived by Getmanov [17] for the f3 > case. In 
addition an expression for the phase 9 which appears in its dual field v, is also obtained 

6 = — Q — arctan ^tan(a) coth ^2y / /jcos(a) ^^ y 2 )) ' 

Respectively for the off-diagonal elements that correspond to the (5 < sector, the set of 
equations is 



0«r-\/\ i\ (S <;■■) = 



(26) 



d^ + ^le^fs+^jil-vv*) = 0, 
that finally produce a different solution 

v = -e in ^cos(a)tanh {^2y[\P\ cos(a) V ^ + zsin(a)^ , (27) 

with Q a real parameter associated with the vacuum of the theory. This is the solution 
that was derived by Lund and Regge fl when considering the (3 < case. 
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A two-soliton solution can be obtained through a non-linear superposition technique. 
Starting from the vacuum of the theory and by the application of the Backhand trans- 
formation twice, a set of parameters {5i,5 2 } is used respectively in each step. The same 
procedure is followed again where the two parameters are used in the opposite order. By 
demanding that the two results are equal, one ends up with an equation for the two-soliton 
solution in matrix form 

92s = CT3 (8l92 ~ $29l) 9vac°3 fal^f 1 ~ S^ 1 ) ■ (28) 



The matrix field variables are of the general form of fllPf) , with elements 

cos(a fc )A^ i; exp(2iv / ^sin(a fc )e fe ) 

u k = 7 \ , (29) 

cosh \ 2y/~/3 cos (afc)Efcj 



Vk = — e in ( cos (ofe) tanh ( 2a/ /3 cos (afc)Efc ) + i sin (a^) ) . (30) 



The identification one must make is: 

e *4(* + £) i+ ^-£> ■ (32) 

where is a total phase. As expected g<i s has the same general form of equation ( [TDD 

92s=( U2s ~ W l s ) • (33) 
The two-soliton solution and its complex conjugate are given by the diagonal elements of 

92s 

u = (-5ml + 5 2 v*) e iQ (jhui - 5 2 u 2 ) + (-to + jggi) e~ in (-61V1 + 5 2 v 2 ) ^ 
5 2 + (-uju 2 - u* 2 ui - v*v 2 - VgVi) 5 2 5 x + 5 2 ' 

u * _ (Siv 2 + 5 2 vi) e- iQ (5ml - hu* 2 ) + (-5m 2 + 5 2 u\) e in (-5ml + 5 2 v* 2 ) 
5 2 + {-u*u 2 - u\u x - v\v 2 - VgVx) 5 2 5 x + 5 2 
while the off diagonal elements represent the dual field and its conjugate 

v _ (Siul + 5 2 up e in (5mi - 5 2 u 2 ) - (-to + favi) e~ in (-5m + 5 2 v 2 ) 
5 2 + (-u*u 2 - u%u x - v*v 2 - v^vi) 5 2 5i + 5 2 

v * = (-5 lU2 + 5 2 m) e~ in (5ml ~ jgn|) ~ (Sm 2 + e in (-5m* + 5 2 v* 2 ) (3J) 

The expressions above represent two-soliton solutions to the equation of motion and are 
related through the duality transformation of (|I5|). 

Multi-soliton solutions can also be obtained by following the same technique. Instead 
of the vacuum solution, one can start from any given n-soliton solution S n and add solitons 
through the method described above, ending up with a S n+2 solution. Nevertheless, multi- 
soliton solutions for this model are quite large and their calculation is beyond the scope of 
this paper. 
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2.4 Soliton - Antisoliton duality 

In this section we argue that soliton-soliton and soliton-antisoliton solutions presented 
in the literature by previous treatments, do not represent distinct classes of solutions. 
Charged solitons are non-topological solutions, therefore a distinction between a soliton 
and an antisoliton is impossible. On the other hand chargeless solutions may be realised 
as topological solitons and identified with the sine-Gordon solitons. 

The sine-Gordon theory appears as the limit of the CSG model when the charge pa- 
rameter a is set to zero. We can substitute in the equation of motion of (|2p 

u = sm0e 2ir > , (38) 

where rj is now a constant to recover the sine-Gordon model in the usual form 

dl4> -dlcf) + 2(5 sin 20 = . (39) 

The sine-Gordon theory has topological solitons (both kinks and antikinks) interpolating 
between its degenerate vacuaa. In contrast the CSG theory has a single vacuum for (3 > 0, 
and therefore its solitons are not topological in nature, but are stable because of integra- 
bility alone. The topological nature is hidden within the mapping of (|38|) and one has to 
be careful when trying to recover the sine-Gordon soliton as a limit of the CSG theory. 
Nonetheless a subtle remnant of the topology survives the mapping to the complex sine- 
Gordon theory. To see this consider how a SG soliton is mapped to CSG soliton. This is 
shown in Fig.|l|. 







Jt/2 


/ * 




\ u 








Figure 1: The and u solitons 

Consider now how the potential term behaves as x increases for the single soliton 
solution. As a function of u we can express the potential as 

sin 20 = 2 sin cos = ±2u\/l — uu* , (40) 
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where r] has been ignored as a total phase. Note that we should take the branch cut with 
opposite signs on each side of the point <fi — |, u = 1. We shall see that the changing sign 
of the branch cut for a chargeless soliton will be important when we come to consider the 
theory with a boundary. 

In some sense the topology of the chargeless w-soliton is embedded in the the branch cut 
that appears at the singular point u — 1. The choice of branch corresponds to a different 
vacuum for <fi and therefore to a different topological charge. 

However, when the charge parameter a is not zero, then the w-soliton does not reach 
the sick point u = 1 and remains non-topological. In this case no real distinction can 
be made between a soliton and an antisoliton. In the sine-Gordon case, the antikink 
solution is derived from the kink by changing the sign of the parameter 5 of the Backhand 
transformation. This effectively corresponds to a parity and time reversal transformation 
which finally produces an antikink solution. Examining ( p9"D we see that in the CSG case 
this change actually leads to the complex conjugate solution, by changing the sign of the 
complex phase. The change of sign in both t and x, can be cancelled by taking the charge 
parameter a — *> —a. It is thus clear that instead of changing the sign of S, one could 
effectively change the sign of a to derive an antisoliton. Since the soliton solution is a 
smooth function in a, the antisoliton is not a distinct object but can be identified with the 
soliton itself. 

This also has an effect on the two-soliton solution. If we follow the same steps as in 
the sine-Gordon two-soliton solution then the solution u 2s of (|34] ) corresponds to both a 
soliton-soliton and a soliton-antisoliton solution depending on the choice of sign for the 
Backlund parameter <5 2 ■ The parameter 5 2 can be chosen in such a way as to describe one 
of the following 



1-V 

5i — — (62) — -w - — — soliton-soliton scattering (41) 



Si = (S2) 1 = y ^ — ^ soliton-antisoliton scattering (42) 

Here we have taken the two solitons to have equal and opposite velocity (In general this is 
not the centre of mass since differently charged solitons have different masses, but it will 
be convenient for our discussion when we introduce a boundary later on). Also for reasons 
of simplicity we will refer to the soliton-soliton solution as u ss ( 62 = — l/S\ ) and to the 
soliton-antisoliton as u sa (62 = l/<5i) . 

However since no topological distinction exists between the soliton and antisoliton sec- 
tor, it is possible to find a transformation of the parameters of the solution which effectively 
acts as a change of sign for the parameter 62- In fact, a set of transformations exists that 
maps u ss to u sa but we restrict ourselves to the simplest cases. 

Before introducing the transformation, we need to introduce arbitrary shifts in x, which 
are crucial not only for this mapping but also later when we consider breathers and soliton 
reflections. The shifts appear in exponentials, so it is more helpful to consider the shifts 
in the following forms 

Ki = exp I 2\I~B cos a* Xl _ \ 
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Ji = exp ^2i^Jpsmai-^=== + mA ; i = l,2. (43) 

The parameters ifj, Ji are directly related with both V and a and correspond to the ar- 
bitrary initial positions in x, in the real (Sj) and imaginary phases (Oj) respectively, that 
appear in the one-soliton solution. Specifically the parameter K represents a translation in 
x, while the J parameter represents a phase shift in the internal U(l) space. For reasons of 
simplicity, we include in the definition of J the total phase iV& = exp(iRk) which appears 
in fl29|) . Henceforth these parameters will be referred as phase shifts, since they are directly 
related to the time-delay effect of the scattering process. 

Now that we have defined the arbitrary phase shifts we start with the soliton-soliton 
solution u ss which comes from the two-soliton solution u 2s when we choose 8 2 = —l/8\. 
We consider the following transformation 

a 2 -> -a 2 . (44) 

Although this is enough to change a single soliton to an antisoliton, this is not the case 
for the two-soliton solution. The phase shifts have also to be fixed in a specific way to 
complete the mapping between u ss and u sa 

Ji -> -Ji (45) 
K 2 -> 1/K 2 . (46) 

This effectively changes the sign of 8 2 in the expression u ss converting one of the solitons to 
an antisoliton. In contrast with the single soliton where the antisoliton can not be properly 
defined, in the two-soliton case there is a point of reference. A distinction between a soliton 
and an antisoliton can only be realised as a specific choice of the relative sign between the 
parameters ai, a 2 and V which does not in any case lead to topologically distinct solutions. 

The same mapping between u ss and u sa can also be achieved by making the following 
transformation 

a 2 — > a 2 + 77 , (47) 

which effectively changes the sign of all trigonometric functions involving the parameter a 2 
sending the solution u ss — > — u sa . This transformation will be used again on a later section 
when we come to consider soliton reflections, to demonstrate exactly the equivalence of the 
two sets of solutions. 



2.5 Breather solutions. 



There are conflicting views in the literature concerning the existence of breathers |17], |18 



The problem arises because the transformation V — > iV which is usually used to generate 
breathers from a two-soliton solution traveling with equal and opposite velocities, does not 
necessarily lead to a solution of the equations of motion. While the technique has been 
widely used before on other models, the fact that the CSG equation involves both u and u*, 
implies that naively analytically continued solutions do not necessarily satisfy the equation 
of motion. 
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So it is not clear, for instance, that all the breather like "solutions" of |14| do satisfy the 
CSG equations of motion. However, since the sine-Gordon is embedded in CSG model by 
taking u to be chargeless, the sine-Gordon breather solutions do satisfy the CSG equations 
of motion. In fact a family of charged, complex breather solutions does exist in CSG 
model. Although it is quite hard to actually check if a general breather solution satisfies 
the equation of motion, a trick can be used instead. We consider the two-soliton solution 
of fl5ip and we demand that this solution is even in V so that is effectively a function of V 2 . 



Now the transformation V — > iV, doesn't change the reality properties of the solution but 
simply introduces an overall minus sign into the arbitrary parameter V 2 , which is irrelevant. 
Making the solution even in V, means that a few restrictions have to be imposed. Firstly, 
the charge parameters have to be taken equal or opposite according whether 5i5 2 is plus 
or minus one respectively. Secondly, some of the arbitrary position parameters, have now 
to be fixed. However, up until now all the arbitrary phase shifts that appeared were either 
complex (Jj) or real (Ki) and there was no distinction between the shifts that originated 
from the space or time part of the phase. However, when constructing a breather solution, 
by analytical continuation of the V parameter a separation between the space and time 
shifts is induced. All shifts that associated with space end up as real parameters, while time 
shifts become imaginary. We can restrict ourselves to shifts only in the x direction. One 
could also consider more general phases which are complex and also depend on time and 
the parameter V. These however correspond to either U(l) rotations or time translations 
which make their use obsolete. The arbitrary shift parameters are now both real 



K s = exp \2J(3 cos a 



VT+v 2 ', 



J s = exp ( 2\fd sin a s . =Hf — | ; s — 1,2 . (4£ 



and should be compared with the general form of (|43|) . In order to make a breather solution 
from the soliton-soliton case the following relations are required 

Ki = ±Tr and J i = t\- , (49) 

where the signs in these equation are correlated. 

It should be noted that more breather solutions may exist. It is possible that through 
certain restrictions a more general breather solution can be obtained, but a direct confir- 
mation through the equations of motion is rather difficult. 



2.6 Collapse of a Breather. 

An analysis of the quantum CSG model |18|] suggests that the soliton can be identified 
with the elementary particle since the vacuum of the theory and the one-soliton are not 
topologically distinct solutions. Evidence of this conjecture exists even in the classical 
picture. From our experience with the sine-Gordon model, we would expect to identify the 
particle with the lowest energy breather solution. It would seem to follow that the breather 
whose energy and charge correspond to that of a single particle should be equivalent to a 
single soliton. This remarkable fact can be shown as follows. 
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We consider the static single-soliton solution 

cos(a) exp(im sm(a)t) 

U s = — -—— r , (50) 

cosn(m cos(a) [x + xq) 

where m = 2yf]3. The mass of the static soliton is given by fllil) which after integration 
gives 

M s = 4mcos(a s ) , (51) 
while the charge of the soliton is given by 

Q s = i jud^i — u^u_ = 4 / ign [ fls ]^ _ a \ ( 52 ) 
J 1 — uu V 2 / 

The mass of a breather is twice the mass of a single soliton solution at velocity V, which 
has been analytically continued V — > iV 

8mcos(a B ) 

M B = - . (53) 

\jl + V B 2 

The breather solution is effectively constructed from two one-soliton solutions, each with 
charge 

Q B = 4 ^sign[a B ]^ - a B ^j ■ (54) 

In order to have a chance of identifying the breather with the soliton, we demand that the 
mass of a breather is equal to the mass of a static single soliton and that their charges also 
coincide 

M S = M B , Qs = 2Q B ■ (55) 
From the above relations, one can solve for the parameter Vb 



Vn 



(56) 
\ cos(a s ) 



If this value is substituted into the breather, then the solution collapses to a static single- 
soliton carrying double the charge Qb- In other words, the single-soliton can always be 
considered as a bound state of two single-solitons carrying half the charge. The argument 
can be used recursively so that a soliton can be regarded as an infinite collection of solitons 
carrying fractions of the original charge. At each level a soliton is identified with a breather 
emerging out of a soliton pair of half the original charge. In the classical picture this process 
can be carried out indefinitely, but in the quantum case the finite character of the mass 
states restricts this procedure. 

This is not surprising since the static single-soliton of fl50D can be viewed as a bound 
state due to the oscillation effect which creates a breather-like behaviour. This is consistent 
with the fact that any breather can collapse to to this solution when the parameter V 
is properly fixed. It can therefore be realised as a breather solution after the collapse, 
exhibiting all of its former properties. 
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One point that has to be emphasized is that breathers constructed with the method 
described in the previous section are not chargeless. This is due to the fact that the 
choice of the charge parameters Oj is such that both solitons that are combined to create 
a breather have the same charge. This is confirmed by the above demonstration in which 
a breather collapses to a single soliton solution which carries double the charge of the 
breather's solitons. Neutral breathers do exist but only at the chargeless limit and can be 
identified with the breathers of the simple sine-Gordon theory. 



3 Reflections. 

In the following sections we consider the effect of introducing a boundary into the CSG 
model. Although we are free to add any boundary potential to the Lagrangian, we choose 
to investigate those potentials and their corresponding boundary conditions which preserve 
integrability. Such choices allow the non-perturbative solution of the model which is one 
of the main motivation for its study. Once such suitable conditions are introduced we 
examine the scattering of particles and solitons, and determine the necessary conditions 
for the existence of boundary bound states. 

3.1 Abelianisation of the Lax pair and conserved currents. 



We shall consider a boundary condition to have preserved the integrability of the CSG 
model, if we can still construct an infinite number of commuting conserved charges. In 
contrast with the theory in the bulk, the introduction of a boundary destroys the translation 
invariance of the model but preserves the time translation invariance. It is thus expected 
that the momentum will not be conserved, whilst the energy will. This situation also 
holds for the higher-spin conserved quantities. All energy-like, parity-even quantities can 
be conserved, unlike their momentum-like, parity-odd partners. Nevertheless, since there 
is an infinite number of conservation laws, the main goal would be to concentrate on the 
conservation of the parity-even quantities. 

The presence of the spectral parameter A in the Lax pair of (|12|) implies the existence 
the infinite conserved currents in the bulk that can be determined through the method 
used by Turok and Olive [2T|. This is achieved by performing a gauge transformation U 



A = UAU' 1 + dUU- 1 , (57) 

in such a way that the commutator of the transformed gauge fields A and A of the Lax 
pair to be zero. The equation of motion becomes 

d (A-A) = d 1 (A + A) , (58) 

where the normal time and space derivatives are used. In the theory in the bulk we 
integrate over x. If Q(X) = J^°°(A — A)dx, then 



d Q(A) = / ~ d (A - A) dx = /"~ di(A + A) dx = \(A + A)} ^ . ( 59 ) 



dt 
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Since at infinity the fields are taken to vanish so that A and A approach a fixed value, 
it follows that the right-hand side vanishes. As A and A can be expanded as an infinite 
Laurent series in A, the coefficients of each power of A, provides us with an infinite number 
of conserved charges. 

When a boundary is introduced the left-hand side involving the spatial derivative does 
not vanish since now the integration takes place over the semi-infinite interval. Instead one 
is left with an equation of the form 

y° d (A -A)dx= [{A + A)] x _ q (60) 

where the left-hand side is evaluated at the boundary. Instead of demanding that the 
right-hand side vanishes, we instead ask that it can be expressed as a total time derivative 
with the help of suitable conditions, thus leading to a conserved quantity. 

We begin by finding explicit expressions for "low-spin" conserved charges of the CSG 
model in the bulk by solving for the abelianizing gauge transformation U order by order 
in the spectral parameter. 

Let U be a general real SU(2) matrix, with det(U) = 1. The diagonal elements of U 
can be taken equal due to residual gauge freedom which leave A and A in an abelian form. 
Thus U takes the form 

^TtM*?)' (61) 

where \ is a function of the fields and should not be associated with the space variable. 
We demand that U diagonalises both A and A at the same time. The transformed fields 
lie both in the 03 direction and the non-zero diagonal elements can be identified with the 
conserved currents. Taking A to be 



A = , (62) 




with A = (3X and E = i{u*dv* — v*du*), we demand that the non-diagonal part of a 
vanishes 



2iA X + X 2 E* + E + d X = , 

2iAx + fE* + E + dx = . 
The conserved quantities can also be written in terms of x an d X 



(63) 



T .,1-XX XE +xE Xdx-Xdx (C , A \ 

J = -*A- — + — — . (64) 

1 - XX 1 - XX 2(1 - XX) 

The same matrix U, should also diagonalise A, which is given by ( pT3[) 

A = ( P* -D ) ' (65) 

where D = uu* — vv* and P = —2iu*v*. The choice of E, P and D is not accidental. 
They actually represent the electric field, the polarization and the population inversion 
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field variables respectively, when this theory is used to describe the propagation of optical 



pulses in a non-linear medium [10|. When U acts on A, we again demand the off diagonal 



parts to vanish. Examining the matrix explicitly yields 

l -(-2D X + P-P*X 2 ) + d X = 



l -(2D X + P* -Px 2 ) + d X = 0. 



(66) 



It is easy to see that these equations are equivalent to equations 

The diagonal part yields the other component of the conserved current 

2X\ I-XX 1 -XX J 2(1 - xx) 

In order to solve the two sets of equations (^) or equivalently (0), we consider an expan- 
sion of x an d X i n powers of A 



= Xl X2 X_3 _ 

X A A 2 A 3 



- _ Xl , X2 , X3 , 

X A A 2 A 3 



(6£ 



The coefficients x« and Xi can be determined by direct substitution into ( p3|) and (J66Q , and 
by demanding that the coefficients in all powers of A vanish. Up to order 0(A -2 ) one finds 

X = (A) E + (^fdE + (^) 3 (E 2 E* + d*E) + ... , 
X=(^)E* + (^) 2 dE* + (^) 3 (EE* 2 + d 2 E*) + ... , 

Now that x and x have been defined, we can also express the conserved quantities as a 
series in A. Each order of A, provides a conserved quantity and since the series of A in 
X and x does not terminate, we thus have an infinite number of conserved quantities as 
expected from the integrability of the CSG model. The two components of the conserved 
current up to 0(X~ 2 ) can be read off as coefficients in the following expansion of J and J 

J = -XP - ^EE- (i) - ^(ESE- - E-OE) (±) + ... , (70) 

and it can be checked that this current is conserved explicitly from the equation of motion. 

In the above we have constructed conserved currents that lead to conserved charges in 
the bulk. However, as we have previously argued, conserved charges on the half line are 
expected to take the form of an integral over a parity- even conserved current. The conserved 
currents above are neither parity even or odd. To rectify this we note that our system of 
equations and constraints possess a Z2 invariance involving parity transformations which 
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can be used to construct a "reflected" set of conserved currents. The "reflected" set of 
conserved currents is easily obtained through the substitution d — > d in the expressions 
(JfO|) including those derivatives involved in the definition of E. The new set of currents J, 
J, can now be combined with the former set to produce pure parity odd and even currents. 

In the presence of a boundary only parity even quantities are conserved. The desired 
form of the equations to emerge is 

d (parity even) = d\ {parity odd). (72) 

By combining the two sets of currents one can separate the odd and even quantities for all 
powers of A. 

d [(J + J) - ( J + J) J = 8 1 [{J - J) + ( J - J)] (73) 
We examine the A -1 term in the expansion which gives 

d [EE* + EE* + 2(3(D + D)) = d x (EE* - EE* + 2/3(D - D)) , (74) 

where E = E(d — > d), etc. After integration over the semi-infinite interval, the right hand 
side representing the parity odd part is 

flbW(»,«0 = G^V) 9 u + a «* • (75) 

V 1 — uu* / V 1 — uu* / 



This is a total derivative provided that 

2d x u* _ dW 2d lU _ dW 
1 — uu* du ' 1 — uu* du* 

The conserved quantity at hand, in terms of u and u*, is then 



(76) 



„ = £ (2 lft ^'_y + - 1)) d, - [W|„„ , (77) 

Since this quantity actually represents the energy of the system, W can be identified with 
the energy contribution of the boundary term. 

When constructing the odd and even quantities of the A -2 term, one ends up with 

do (~(E*dE - EdE* + EdE* - E*dE) - f3(E*P - EP* + E*P - EP* 
dx (^-{-E*dE + EdE* + EdE* - E*dE) - f3{E*P - EP* - E*P + EP*)J (78) 

Once more the parity-odd right hand side which after integration yields 

^doudodiu* ^d u*d diu ^{diudiu* + doudou*){udiu* — u*diu) 
1 — uu* 1 — uu* (1 — uu*) 2 

-4»-«-^ + 4^-4^ , (79) 

1 — uu 1 — uu 
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should be written as a total time derivative in order to force the currents to be conserved 
at the boundary. A set of boundary conditions have to be introduced to ensure that his 
is the case (App. |A.1| ). Using the equations of motion, the parity-odd part of (|79|) can be 
written as a total derivative if the following restrictions are enforced 

d\u = —Cu\J\ — uu* , 

(80) 

diu* = —Cu*y/1 — uu* . 

The boundary constant C is defined by the theory and is responsible for the way fields 
react to the boundary. Consistency of the two equations in (|80|) implies that C should be 
considered a real parameter. When one makes the transformation described in (^j, the 
new boundary conditions for the fields and rj are 

d l( f) = -Csin(0) , 8^ = , (81) 

which clearly shows that C has to be real. 

It has to be pointed out that (|StJ) is not the only set of boundary conditions that can be 
derived. A number of isolated "Dirichlet"-like conditions also exists. However, we restrict 
ourselves only to cases where the space derivatives of the fields appear. If we take the field 
u to be real, the system is reduced to the sine-Gordon equation with a boundary condition 
d\(f) = —C sin0 . This is the subset of integrable boundary conditions of the sine-Gordon 
model presenting the Z 2 symmetry <p — > —cf). The corresponding conserved quantity, is 
rather large and is omitted. 

4 Soliton scattering and boundary bound states 

Since the necessary conditions for the integrability of the model have been established, we 
study the scattering of particles and solitons off the boundary. We begin this section with 
the effects of introducing a boundary potential to the vacuum of the theory. We continue 
with the scattering of particles and solitons and derive the phase shifts induced by the 
process. Finally, we investigate the necessary conditions for the existence of boundary 
bound states. 

4.1 Vacuum 

When a boundary term is introduced, the vacuum of the theory that we discussed in sec- 
tion (2.2), does not necessarily remain unchanged. It is exactly this contribution that 
needs to be carefully examined before any statements are made about the minimum en- 
ergy configuration. Although, the vacuum solution of the theory in the bulk is a strong 
candidate, soliton solutions could also be considered in the attempt to both minimize the 
energy functional and satisfy the boundary conditions of (|80|) . 

We begin by first determining the energy contribution of the boundary term. The full 
Lagrangian of the model is now 

C tot = C + C B . (82) 
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The boundary term C B , can be determined by the variation principle of the total action. 
The variation of the £ term yields 



When the Euler-Lagrange equations are used, two terms survive since the model is consid- 
ered in the semi-infinite interval where the fields do not vanish at the boundary 

From the variation of the boundary term one has 

5C B = d -^5u + f^5u* . (85) 
du ou* 

The variation of the action vanishes when the remaining terms evaluated at the boundary 
are forced to cancel. The two interrelated equations that emerge are 

(86) 
(87) 

By substituting the boundary conditions of fl8"0D, these can easily be solved for the 
boundary term 

C B = 2CV1 - uu* . (88) 
We now consider the total energy of the system , now comprising of two parts 

7~Ltot = 7~Lbuik + 7~~(-b , (89) 

where the term Tibuiki represents the energy in the bulk and the second term Tis represents 
the energy contribution from the boundary 



dC 


—diu 


dC B 


d(diu*) 


1 — uu* 


du* 


dC 


-d lU * 


dC B 


d{d\u*) 


1 — uu* 


du* 



U B = -2CV1 - uu* , (90) 

which is evaluated at x — 0. This energy contribution makes the determination of the 
vacuum difficult. The sign of the boundary constant C is not set, which could provide 
either a positive or negative contribution to the total energy of the system. This clearly 
shows that although the original choice for a vacuum should not be discarded, one should 
also consider other static solutions which in conjunction with the sign of C could provide 
a lower energy vacuum than before. 

Apart from the original choice for a vacuum, one can consider static multi-soliton so- 
lutions. We restrict ourselves to one-soliton solutions since experience with similar models 
usually makes multi-soliton solutions unsuitable candidates. 
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When considering one-soliton solutions, one has the equations of the Backlund trans- 
formation (|23|) which are always true to simplify expressions. In particular we first consider 
the TCbuik term representing the energy in the bulk 

H 6u zfc = / <ir h m uu ) , (91) 

J y l — mm* y 

with m = 2y^. By direct substitution of the Backlund equations a simplified expression 
of the bulk energy is acquired 

Hbuik = J dx{2m 2 uu*) . (92) 

When the above expression is integrated throughout space, the result can be identified 
with the mass of the soliton solution u. However, now the integration is over the half line 
and specifically over the [— oo, 0] region. 

The same equations can be used to express the Hb term of (|90| ) . Specifically, the 
boundary constant C is determined by direct comparison of the Backlund equations of 
(^3|) and the boundary condition which appears in (jS0[) 



C=-(s e« e+Q 1 + - c e~^ +n ) ) . (93) 



with 6 given by (|25|). At x = and assuming that V = 0, the above expression simplifies 
to 

C = ± . (94) 

y 1 + tan 2 (a) coth 2 (mcos(a)xo) 

This implies that \C\ < \m\. Since both m and C are defined by the theory, the above 
relation is true only for specific choices of the boundary parameter C. Alternatively, one 
can think of this restriction emerging from the fact that for \C\ > \m\, no choice of xq 
satisfies the boundary condition. 

In the case where we choose u = as a possible vacuum, the only remaining term in 
the total energy is 

Htot = -2C . (95) 

Alternatively, one can consider a one-soliton solution where V is set to zero which appears 
in (0). In this case, both terms of fl89D depend on the initial position of the soliton. 
However, after some calculations, the xq dependence drops out and the total energy is 
given by a the following expression 

'Htot = 2mcos(a) . (96) 

It is far from obvious, which vacuum choice provides the minimum energy configuration. 
To determine this, one has to look at the expression of the boundary constant C in ([94]) . 
This can be rewritten in the following form 

y 2 + ^f^ = cos 2 (a), (97) 



19 



where 

C 

F 2 = sinh 2 (m cos(a)a;o) , y — — • (98) 

m 

In the above relation m and C should be treated as fixed parameters, while F can be varied 
through xo- The left hand side of (|97|) is monotonically decreasing as F increases since 
< y 2 < 1. We observe the following 

cos 2 (a) = 1 when F — > , 

(99) 

cos 2 (a) = y 2 when F — > ±oo . 

This shows that moving the soliton away from the boundary decreases the energy of the 
system. On the extreme case where the soliton is placed at infinity, the model behaves as if 
no soliton exists, and the only contribution is the boundary term which coincides with the 
vacuum solution of u = 0. On the contrary as the soliton is placed closer to the boundary 
the energy increases. The maximum energy occurs when F = at which point cos 2 (a) = 1 
so that H tot = 2m which is greater than the energy H tot = 1C of the u = vacuum. 

Although the choice of vacuum in the bulk seems to be the most suitable choice in the 
boundary case too, one cannot rule out multi-soliton solutions that might provide lower 
energy configurations. This demands tedious calculations and remains as one of the open 
questions for this model. 

4.2 Soliton reflections 

In this section we investigate the reflection of solitons from the boundary. Mathematically 
this can be represented by a two-soliton solution satisfying the boundary condition. One 
of the solitons represents the incoming soliton whilst the other represents the reflected 
one. The point where the two solitons actually meet along the whole line as well as the 
phase shift due to their collision create an overall time-delay effect which can be calculated 
directly through the parameters of the scattering. This time-delay can be attributed to 
the interaction of the soliton with the boundary. 

However, the most difficult step is to determine the restrictions that have to be imposed 
so that the two-soliton solution satisfies the boundary condition 

d\U2s = -Cu 2s \]l - u 2s u* 2s . (100) 

Energy and charge conservation laws demand that both the mass and the charge of the 
soliton are conserved by the boundary. This restricts the choice of the charge parameters 
ai, a 2 to be either equal or opposite. 

Due to the large expressions involved in the calculation, one is forced to expand both 
sides of the equation ( |100| ) to a Taylor series in exponentials of t, and match each term of 
the same order. Each term provides us with an equation involving the boundary parameter 
C . As mentioned in the previous section, the boundary constant has to be a real parameter. 
The real and imaginary parts of the equation yield two constrains on the parameters. 

Let us consider this in more detail. We begin with a two-soliton solution, where the 
parameters are chosen in such a way so as to describe a soliton-soliton scattering. In 
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this case, the charge parameters are taken to be opposite a± = 
parameters to be 5\ = — 1/5 2 - 

Furthermore, we adopt the following parametrisation which is more natural 



ct2 and the Backhand 



Ki 
K 2 



■h 
■h 



V = tanh(tf) 



(101) 



After both sides of the boundary equation are expanded as a Taylor series in time, we can 
discard the imaginary parts from all terms by using the following relation 

sin(C) = - tanh(tf) tan(a) sinh(A). (102) 

When the above equation is used the infinite set of equations collapse to a single constraint 

'cos(C) + cosh(A) \ 



C = m cos(a) cosh($) 



sinh(A) 



(103) 



When the shift parameters are fixed according to the above relations, the two-soliton 
solution satisfies the boundary condition and this process describes a soliton being reflected 
by the boundary. 

The fact that only relative shifts in both normal and internal U(l) space are important 
should be expected from time translational and U(l) invariance of the model. The non- 
topological solitons in the CSG theory are reflected as solitons carrying the same charge 
Q. This is because the boundary potential does not break the U(l) symmetry since it 
depends only on \u\. 



. x 





Soliton-Soliton sector 

Soliton-Antisoliton sector 



Figure 2: C in terms of A 

We shall consider a and $ as fixed parameters and use equations ( |102[ ) and ( |103j ) to 
determine the parameter A in terms of C. By eliminating the parameter (, one ends up 
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with the following quadratic equation for C 



C 2 - 



Cm cos(a) cosh(-#) cosh(A) 



+ (sinh(^) 2 + cos(a) 2 ) m 2 = . (104) 



sinh(A) 



The solutions of the above equation can be plotted to present the dependence on A. 
The plot involves two branches (Fig. |2|) due to the sign ambiguity, which are mutually 
exclusive. The plot shows that a soliton can always be reflected by the boundary. The 
branches meet at the points 



In the limit a — ► the two branches of the plot can be identified with the soliton-soliton 
and soliton-antisoliton sector of the reflection process at the sine-Gordon limit (Fig. |3|). 
For fixed values of d and a = 0, it is the value of the boundary constant C which determines 
whether a soliton is reflected as a soliton or an antisoliton. For C small, a soliton is reflected 
as an antisoliton (Neumann boundary conditions for C = 0), while for C large a soliton 
is reflected as a soliton (Dirichlet boundary conditions for C = oo). For C = mcosh(^) 
the branches do not meet as in the CSG case. This specific value of C corresponds to 
a logarithmic divergence that appears in the classic time delay for the sine-Gordon case. 
These results coincide with the results derived by previous treatments of the boundary 
sine-Gordon model p. 




coth(A) 




(105) 



c 




X 








Soliton-Soliton sector 
Soliton-Antisoliton sector 



Figure 3: C in terms of A for chargeless case 
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4.3 The classical time delay 



The time delay which appears at the scattering of a soliton off the boundary, can be 
calculated directly from the asymptotic values of the solution at t — ±00. We begin with 
the two-soliton solution and change to a frame of reference which moves with the incoming 



soliton (i.e. 



x 



Vt + x. In the limit t 



-00 the solution becomes 



where 



and 



S- 



P sin (a) 



cos a 



lim u 2s 



cos 



(a) e ^ +B ^ 



cosh(P(— x + x\) 



(1 - V 2 )t -Vx-yi , tan(5i 



V sin(a) 



cos a 



P 



m cos(a) 



V 2 



iln — 

2 cos 2 (a) + V 2 sin 2 (a) 



(106) 



(107) 



108) 



v/r^V 2 ' 

The parameters Xi and yi represent regular shifts that were introduced in ([431) . The so- 
lution, as expected, describes a single incoming soliton at early time far away from the 
boundary. 

We repeat the same calculation, but now we change to the frame of reference of the 
outgoing soliton (i.e. x = —Vt) and calculate the limit of the two-soliton solution at 
t = +00 which yields 



where 



and 



A, 



S 4 



P sin(a) 



cos a 



cos(a) e i{M+B2+ ^ 

lim U2s = 7—, — 7 — r r 

t=+co cosh(P(— x — X2) + q) 



[1-V 2 )t + Vx~y 2 , tan(B 



V sin(a) 



cos a 



q 



-In 

2 cos^ l a 



V 2 



2f„\ +V 2 sin 2( a ) 



(109) 



110) 



111^ 



Once again this is a single soliton solution representing the reflected soliton far away from 
the boundary wall. 

The asymptotic solutions S+, S- contain all the information needed to calculate the 
time-delay. The latter is a combination of two separate events. Firstly, a phase shift is 
induced during the scattering of the two solitons. Before the two solitons re-emerge as 
two separate entities, the reconfiguration of the solution creates a phase shift which is 
equivalent to a time delay. Secondly, the centre of mass of the two-soliton solution does 
not necessarily lies at the boundary. This implies that the two solitons actually meet at a 
different point than x — 0. This creates again a time delay which may be either positive 
or negative corresponding to an attractive or repulsive boundary potential respectively. 

Ignoring any interaction between the two solitons, one can project the trajectories of 
S + and S- on x — t diagram and find the point where these cross (Fig. ^). The distance of 
this point from the boundary is proportional to the time delay At which in the diagram 
is given by the distance (AC). The time delay corresponds to the time interval in which 



23 




Figure 4: Scattering diagram 



the soliton appears to be absorbed by the boundary before it reemerges as a well defined 
entity . The two solitons move across the following lines 

S_: t = - (P(-x + xi) + r) , 
S+ : t = ^ (P(x + x 2 ) - q) , 



as dictated by ( p.U9| ) and ( |106| ). The lines cross at 



If r + q\ xi - x 2 r 
x o=2{ Xl ~ x 2 + _ ^J = — 2 P ' ^ ' 



since r = q. The time delay is finally 



Atcsg = — = 77 + Z^7T373 ln I ZZiTTYT T/2 -2/ ^ 1 • ( 113 ) 



V V mV cos(a) \ cos 2 (a) + V 2 si 



sin a 



In the expression above, the first term of the right-hand side represents the time delay 
caused by the non-symmetric character of the solution with respect to the boundary In 
the special case where x% = x%, the centre of mass lies on the boundary and the term 
vanishes. The second term is independent of the initial position of the two solitons or the 
boundary potential and is caused by the phase shift of the scattering process. 

The relative position of the two solitons are however fixed according to the constraint 
equations (|102|) and (103) which ensure that solution satisfies the boundary condition. 



Specifically the parameter A corresponds exactly to the x\ — x 2 difference up to the overall 
factor P. It is thus possible to express the time delay in terms of the boundary constant, 
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by solving the constraint equations and substituting the relative position of the solitons. 
We choose to express the velocity parameter V in terms of the rapidity $ for simplicity 
reasons. After a few straightforward calculations we recover the following expression for 
the time delay 

hiQ x 

&TCSG = ~ ? !/■ i / qx 114 

2m cos^a) smn(#J 

where 

sinh 4 ^ ((cos 2 (a) + sinh 2 (^))m 2 + 2Cmcos(a) cosh(^) + C 2 ) 
(cos(a) 2 + sinh 2 ^)) 2 ((cos 2 (a) + sinh 2 (^))m 2 — 2Cmcos(a) cosh(^) + C 2 

In the special limit of a = 0, the time delay for the sine-Gordon model is recovered 



1 , / , 9 / m cosh(^) + C 

Atsg = ■ uV\ ln tanh W h o n 

msmn(w) y m coshf^) — C i 



(116) 



This is exactly the time delay calculated for the sine-Gordon theory in the presence of a 
boundary for the restricted class of boundary conditions which admit = as a vacuum 
to which the chargeless limit of CSG correspond. 



4.4 Boundary bound states 

In this section we examine the spectrum of bound states. Once again, for the the boundary 
condition to be satisfied we need to restrict some of the parameters in the solution. 

The simplest bound state that we can have is the static single soliton that was intro- 
duced in (p0|) . The solution is not really static, as the imaginary phase survives the setting 
of the speed parameter V to zero. The solution is static only in the sense that the centre of 
mass doesn't translate in the x direction, although the wave oscillates with fixed angular 
velocity u = msin(a). 

When a boundary is introduced a static soliton can satisfy the boundary condition for 
|C| < |m| when its position is fixed according to equation (p4[). At the chargeless limit 
any time dependence vanishes and the solution collapses to a static single soliton of the 
sine-Gordon theory, fixed at the boundary. 

Breathers that have been constructed by the method described in section ( |2.5|) can 
also be shown to satisfy the boundary condition. The condition that C is real still holds. 
However, all the arbitrary phase shifts are now real numbers and constrained. We ex- 
amine breather solutions that emerge from the soliton-soliton case. Just as before, the 
solution does satisfy the boundary condition with some restrictions involving the arbitrary 
parameters. Once more a Taylor expansion of the boundary equation is needed. The 
parametrization used in this case is 

K x = e x , Jt = e c , V = tan(tf) , (117) 

while the parameters K2 and J2 have been been properly fixed so that this is a breather 
solution. The first restriction needed for the solution to satisfy the boundary condition is 

sinh(2C) = - tan(a) tan(tf) sinh(2A) . (118) 
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The parameter i9 plays the role of the rapidity, which has now been analytically continued. 
The second restriction which completes the necessary requirements for a boundary bound 
state is 

r n ^ f cosh(2C)-cosh(2A) ^ 

C = mcos(a) cos($) I sinh ( 2 A) J' ^ ^ 

Both relations can be recovered by analytical continuation of the corresponding relations 



of ( |102| ) and ( |103| ) after the necessary restrictions for a breather solution have been already 
taken into account. It is instructive to examine the relation between the parameters C and 




Figure 5: C in terms of A for Breather case 

A (Fig. |5|) since it provides valuable insight to the structure of bound states. There are 
two distinct regions of values of C that do not correspond to any bound state. This regions 
are defined by the limit values 

C = ±m (cos(a) cosh(fl) ± sin(a) sinh(0)) . (120) 

At the chargeless limit, the regions collapse to the single values 

C = imcoshfl , (121) 

which coincided with the logarithmic divergence appearing in the time delay for the soliton 
reflection. 

4.5 Particle Reflections. 

In this section we consider the spectrum of particles and their reflection factors in the 
presence of a boundary. 
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For small fluctuations around the vacuum u = 0, the boundary condition ( j30|) becomes 



die(x,t) = -Ce(x,t). (122) 

In order to calculate the reflection factor when particles bounce of the boundary wall, we 
substitute in the last relation the particle solutions presented in (|19|). The constant A of the 
right propagating waves is taken to be one, since it has to do with the characteristics of the 
particle beam. The reflection factor is identified with the constant B, which corresponds 
to a phase change as the particles encounter the boundary 

B=p%. (123) 
ik — C 

The reflection factor, as expected, depends on C which as stated before, appears as a 
free, real parameter in the boundary condition. For C = 0, the reflection factor is equal to 
B = 1 and no phase appears between the two waves upon their scattering off the boundary. 
This is consistent with the fact that the boundary term is proportional to the boundary 
constant, so when C is set to zero, the boundary term vanishes. 

Particle solutions can be related to bound states through the pole appearing in B. 
Indeed one may choose k = —iC and apply this to a solution of the form ^e(x,t). The 
remaining terms depend explicitly on the boundary constant 

e(x,t) = e -i{^ m2 ~ CH +~ iCx ) . (124) 

This solution is square integrable only for a specific range of values for the boundary 
constant. Specifically if C is positive then the solution is not square integrable since it is 
exponentially increasing as x —>■ oo. 

When — m < C < , then e(x, t) represents a square integrable exponentially decreasing 
solution as x — > oo. It oscillates with constant angular velocity to = \J m 2 — C 2 and is 
therefore a stable bound state. It can also be viewed as the tail of a static one-soliton 
solution satisfying the boundary condition, with the parameters adjusted in such a way its 
centre of mass goes to positive infinity. Examining the condition (|94]) for the static soliton 
to obey the boundary condition, this limit can be achieved as Xq — > oo, i.e. we must take 
the charge is such a way that C = —m cos a. 

Finally in the region C < —m, the solution can increase exponentially in time. This 
shows that the vacuum solution u = is no longer stable. In fact the particle behaviour 
which corresponds to a small perturbation around the vacuum seems to be ever increasing. 
This instability can be understood through a rather impressive mechanism in which a 
chargeless soliton is emitted from the boundary, effectively changing the value of C so that 
u = is now stable. 

Recall from section (|2.4j) that for a chargeless soliton we should take the opposite sign 



for y/1 — uu* on each side of the centre of the soliton where \u\ = 1. The instability can be 
viewed as a left moving chargeless soliton which approaches the boundary from x = oo. In 
the beginning while the soliton is far away from the boundary u = so that the boundary 
potential of is H B = —IC . As the centre of the soliton passes through x — 0, the sign 
of the square root in the boundary potential changes. As the soliton moves to x — — oo, 
u returns to near the boundary but now we take the boundary energy with the opposite 
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sign Hb = 2C. Effectively the sign of C has been flipped to a positive value. The energy 
released from the boundary is AC > 4m, which is greater than the rest mass of a single 
chargeless soliton. At C = —m, the soliton is emitted with infinitesimal velocity. As C 
decreases, more energy is given up by the boundary and the soliton can be emitted with 
larger V. This process agrees with the infinite time-delay effect which was encountered in 
the soliton reflections section. The soliton emission represents the time reversal picture of 
that effect in the chargeless soliton case (Fig. [3]). It follows that we need never consider 
the situation where C < — m. 



5 Discussion 

The CSG model is one of the simplest generalisations of the sine-Gordon theory, but 
nonetheless has a rich and fascinating mathematical structure. 

In the first part of the paper we have examined the spectrum of the theory in the 
bulk and written down explicit two-soliton solutions within the framework of the matrix 
potential. We also demonstrated how to construct breather solutions in an elegant way 
avoiding the problems that arise by the analytical continuation of the parameter V. There 
are two ways in which soliton solutions of the CSG model differ qualitatively from those of 
the Sine-Gordon model. Firstly, the solitons can be charged, so that whilst there centre of 
mass remains fixed, the solutions are not static. The second feature is that the CSG does 
not possess degenerate vacua, and so the solitons are non-topological. There is therefore 
no distinction between solitons and antisolitons which can be interchanged by a continuous 
variation of the charge parameter a. Nevertheless, the topological nature of the sine-Gordon 



theory can be recovered as the choice of branch cut of Jl — \u\ 2 in the chargeless limit 
as was demonstrated in section ETJ] • A direct consequence of the non-topological nature 



of the CSG soliton is that the breather solution can collapse to a single soliton when the 



parameter V is properly fixed, tying in with the picture presented in [113] that particle and 
solitons can be identified in the quantum limit. 

In the second part of the paper we introduced a boundary term in the CSG Lagrangian 
and demanded that the system remains integrable. First we constructed low-spin conserved 
quantities of the theory using abelianisation of the Lax pair, and then derived suitable 
boundary conditions in order to preserve these. In the presence of a boundary, we examined 
the vacuum structure and showed that the bulk vacuum u = remained the true vacuum 
in the boundary case. 

Soliton reflections off the boundary were also studied and the necessary constraint 
equations were written down in terms of the phase shift parameters. The set of equations 
was derived by demanding that the two-soliton solution satisfies the boundary condition. 
Moreover the time delay induced by the scattering process was calculated in terms of the 
boundary constant C and was found to coincide in the chargeless limit with the time delay 
of the sine-Gordon theory. 

Finally we looked for classical solutions corresponding to boundary bound states. We 
found that it was possible to construct both bound soliton and bound breather solutions. 
We also found a bound state in the particle spectrum. This was unstable when parameter C 
associated with the boundary energy was in the range C < —m. In this case the boundary 
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emits a chargeless soliton, effectively changing the sign of the boundary parameter C to 
C > m. 

We end our discussion by pointing out a few aspects of the model that appear quite 
interesting and deserve further study. As mentioned in the introduction the CSG theory 
was used to generalise the existing field theory approach of optical pulse propagating in a 
non-linear medium. A physical interpretation of the results appearing in this paper would 
be extremely interesting, especially the physical meaning of breather solutions and their 
application to physical geometries. 

An obvious extension of our results is to consider the quantum case of the boundary 
CSG model. The S-matrix for the model in the bulk, which corresponds to perturbed Z n 
parafermions, is known. It would be interesting to see if one can find a quantum reflection 
matrix, compatible with this S-matrix and with the classical results presented in this paper. 
As the simplest case in the family of homogeneous sine-Gordon theories, the results might 
shed light on the more complicated models in the family. 
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A Appendix 

A.l The boundary condition from the A 2 term 

In order to construct an infinite number of conserved quantities when a boundary is intro- 
duced, we need to express the parity odd terms of the equations of motion as total time 
derivatives. The corresponding term of order A 2 is in terms of the fields u and u* 

o (d x u*) (dlu + dfu) A (d u*) d x d u ^ (diu) (dlu* + dfu*) ^ (d u) d^u* 



1 — uu* 1 — uu* 1 — uu* 1 — uu* 

u* (diu) 2 diu* u (diu*) 2 diu u (d u*) (d±u*) d u u* (d u) 2 diu* 



UU*) (1 — UU*) (1 — UU*) (1 — UU* 

w(<9 u*) 2 <9iw u* (d u) (d x u) d u* # 

-2-^ ^r + 4 n ^ +±(ud 1 u -u d x u)P . 

(1 — uu*) (1 — uu*) 

Since we need the above expression to be a total time derivative we can eliminate any 
second order spatial derivatives of the fields by using the equations of motion of (||) 

^doudodxu* ^d u*d diu ^(diudiu* + d ud u*) (udiu* - u*diu) 
1 — uu* 1 — uu* (1 — uu*) 2 

-A^ud^-u^+A^^-A^^ , (125) 

1 — uu* 1 — uu* 

We take advantage of the fact that we are free to add total time derivatives on this expres- 
sion, since this represents a conserved quantity. The expression simplifies significantly by 
adding the following term 

diu*d u — diud u 



V 1 — uu* I 



which yields 



1 (<9i u*) dlu | i (di u) d 2 u* | 4 (d u) d x d u* ^ (d u*)did u 
1 — uu* 1 — uu* 1 — uu* 1 — uu* 

^ (-udi u* + u* di u) ((di u*)d 1 u+ (d u) d u*) 

(1 — uu*) 2 

+4 (-vdx u* + u*d l u){3 . (127) 
We are looking for boundary conditions that are of the form 

d 1 u = F(u,u*) , d lU *=G(u,u*) , (128) 

where F and G are functions of the fields not involving derivatives. By direct substitution 
of the above into ( |127| ) we get 

a (2j|Ct - uu*) + Gu*) (d u) 2 (d u*) 2 (2^(1- uu*) + Fu 



uu*) 2 (1 — uu*) 2 

Jd u*)(-^ + §^)d u a (-uG + u*F)GF 

o h 4 ^ 

1-ttu* (1-uu*) 

+4 (-uG + u* F) (3 . 
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The expression above does represent a total derivative when all terms are forced to vanish 
by selecting suitable functions F and G. The two separate differential equations that 
appear involving the undefined functions 

dF 

2(1 -uu*) — +uF = , 
ou* 



BC 

2(1 -uu*)— + u*G = , 
ou 



can easily be solved to yield 



F(u,u*) = S^y/l-uu* , G(u, u*) = S 2 {u*Wl-uu* . (129) 
In addition, the last two terms in (|129|) imply that 

F = -G , (130) 



■IT 



Using the above relation and solutions of (|129| ) into the remaining terms of ( |129| ), we can 
determine the remaining undefined functions Si and SV The final form of the boundary 
conditions are 



diu = —Cuy/1 — uu* , 

(131^ 



d\u* = —Cm* yl — uu* 

where C is a real constant. 



References 

[1] P. Bowcock, E. Corrigan, P. E. Dorey, and R. H. Rietdijk. Classically integrable 
boundary conditions for afline toda field theories. Nucl. Phys., B445:469-500, 1995. 

[2] H. Saleur, S. Skorik, and N. P. Warner. The boundary sine-gordon theory: Classical 
and semiclassical analysis. Nucl. Phys., B441:421-436, 1995. 

[3] E Corrigan. On duality and reflection factors for the sinh-gordon model with a bound- 
ary. Int. I. Mod. Phys., A13:2709-2722, 1998. 

[4] Fernando Lund and Tullio Regge. Unified approach to strings and vortices with soliton 
solutions. Phys. Rev., D14:1524, 1976. 

[5] K. Pohlmeyer. Integrable hamiltonian systems and interactions through quadratic 
constraints. Commun. Math. Phys., 46:207, 1976. 

[6] Fernando Lund. Example of a relativistic, completely integrable, hamiltonian system. 
Phys. Rev. Lett., 38:1175, 1977. 



31 



[7] Carlos R. Fernandez-Pousa, Manuel V. Gallas, Timothy J. Hollowood, and J. Luis 
Miramontes. The symmetric space and homogeneous sine-gordon theories. Nucl. 
Phys., B484:609-630, 1997. 

[8] Carlos R. Fernandez-Pousa and J. Luis Miramontes. Semi-classical spectrum of the 
homogeneous sine-gordon theories. Nucl. Phys., B518:745-769, 1998. 

[9] H. J. de Vega, J. Ramirez Mittelbrunn, M. Ramon Medrano, and N. Sanchez. The 
general solution of the 2-d sigma model stringy black hole and the massless complex 
sine-gordon model. Phys. Lett, B323:133-138, 1994. 

[10] Q-Han Park and H. J. Shin. Field theory for coherent optical pulse propagation. Phys. 
Rev., A57:4621, 1998. 

[11] S.L. McCall and E.L. Hahn. Self-induced transparency by pulsed coherent light. Phys. 
Rev. Lett, 18:908, 1967. 

[12] Ioannis Bakas. Conservation laws and geometry of perturbed coset models. Int. J. 
Mod. Phys., A9:3443-3472, 1994. 

[13] Q-Han Park and H. J. Shin. Classical matrix sine-gordon theory. Nucl. Phys., 
B458:327-354, 1996. 

[14] Q-Han Park and H. J. Shin. Duality in complex sine-Gordon theory. Phys. Lett., 
B359:125-132, 1995. 

[15] V. G. Knizhnik and A. B. Zamolodchikov. Current algebra and Wess-Zumino model 
in two dimensions. Nucl. Phys., B247:83-103, 1984. 

[16] H. J. de Vega and J. M. Maillet. Semiclassical quantization of the complex sine-gordon 
field theory. Phys. Rev., D28:1441, 1983. 

[17] B. S. Getmanov. New lorentz-invariant system with exact multisoliton solutions. Sov. 
Phys. JETP Lett, 25:119, 1977. 

[18] Nicholas Dorey and Timothy J. Hollowood. Quantum scattering of charged solitons 
in the complex sine- gordon model. Nucl. Phys., B440:215-236, 1995. 

[19] Q. H. Park, Phys. Lett. B 328 (1994) 329 ||arXiv:hep-th/9402038|| . 

[20] S. Ghoshal and A. B. Zamolodchikov, Int. J. Mod. Phys. A 9 (1994) 3841 [Erratum- 
ibid. A 9 (1994) 4353] ||arXiv:hep-th/9306002|| . 

[21] D. Olive and N. Turok. Local conserved densities and zero curvature conditions for 
toda lattice field theories. Nucl. Phys., B257:277, 1985. 



32 



